On the magnetic ground state of an experimental S — 1/2 kagome antiferromagnet. 
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We have carried out neutron powder-diffraction measurements on zinc paratacamite 
Zn :E Cu4- 2 ,(OH)6Ci2 with x = 1, and studied the heat capacity in fieids of up to 9 T for 0.5 < x < 1. 
The x — 1 phase has recently been shown to be an outstanding realisation of the S = 1/2 kagome 
antiferromagnet. A weak mixing of Cu 2+ /Zn 2+ between the Cu and the Zn sites, corresponding to 
~ 9% of all Cu 2+ for x = 1, is observed using neutron diffraction. This "antisite disorder" provides 
a consistent explanation of the field dependence of the heat capacity for 0.8 < x < 1. From compar- 
ison of the derived Cu 2+ occupancy of the Zn sites for x = 0.8 ... 1 with the magnetic susceptibility, 
we argue that for x = 0.8 ... 1 zinc paratacamite is a spin liquid without a spin gap. The presence 
of unpaired but nevertheless strongly interacting spins gives rise to a macroscopically degenerate 
ground state manifold, with increasingly glassy dynamics as x is lowered. 

PACS numbers: 75.40.Cx,75.45.+j,75.30.Hx 



Physical realisations of the 5=1/2 kagome Heisen- 
berg antiferromagnet have been long sought after be- 
cause it is expected that the ground state of this sys- 
tem can retain the full symmetry of the underlying ef- 
fective magnetic Hamiltonian [l|, 0] ; the geometry of the 
kagome lattice frustrates the classical Neel antiferromag- 
netic ordering, and no symmetry-breaking transition is 
expected even at T = [3|, 0, H, IS • It has been sug- 
gested that even in the thermodynamic limit the symmet- 
ric quantum-mechanical electronic ground state is pro- 
tected from quantum-mechanical dissipation § by a gap 
between the non-magnetic ground state and the lowest 
magnetic (triplet) excitations 0, In practice, how- 
ever, a gap of order J/20 @ may be too small to expect 
a non-magnetic ground state in real materials. 

Shores et al. llj have shown that in the Cu salt her- 
bertsmithitc 

P3 (ZnCu 3 (OH) 6 Cl 2 , depicted in the inset 
of figure [1]) antiferromagnetically coupled Cu 2+ ions are 
located at the vertices of a kagome lattice. Muon experi- 
ments have shown that the ground state of this system is 
either paramagnetic or (quantum) spin-liquid. Almost no 
muon relaxation was observed even at 50 mK |l3| , despite 
the large Weiss temperature W w -300 K jlll.ll4|. Sep- 
arating the kagome layers are Zn sites of Oh symmetry, 
which can also host Cu 2+ ions to form the zinc parata- 
camite family of stoichiometry Zn x Cu4_ x (OH)6Cl2 with 
0.3 < x < 1. For Zn 2+ stoichiometries x < 0.3, the Zn 
site is mainly occupied by Jahn- Teller active Cu 2+ ions 
and becomes angle distorted. At this point the symme- 
try of the lattice changes from rhombohedral (x > 0.3) to 
monoclinic, forming the end-member clinoatacamite [la ]. 
Due to the strong Jahn- Teller distortion of the Cu sites on 
the kagome lattice, the mixing of Cu 2+ and Zn 2+ between 
the Cu and the Zn sites (antisite disorder) in the x = 1 
phase can be expected to be low, but exactly how low 



has not previously been measured. From the magnetic 
susceptibility of samples with x < 1 it is clear that the 
Cu 2+ ions on the inter-plane Zn site are only weakly cou- 
pled to the kagome layers, and it has been suggested that 
the divergence of the magnetic susceptibility for x = 1 at 
low temperatures can be explained by antisite permuta- 
tions of 6-7% of the Cu 2+ ions with - 19% of the Zn 2+ 
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FIG. 1: Section of the neutron diffraction pattern of 
ZnCu3(OD)6Cl2 at 10 K, and the result of the structure re- 
finement (residuals x 2 = 13.95, R p = 2.83%, R wp = 2.73%). 
The inset shows two kagome layers with Cu sites (black), and 
the interplane Zn site (grey) of the zinc paratacamite struc- 
ture. The antisite disorder means that a fraction of the Cu 2+ 
ions interchange with (for x = 1 an equal number of) Zn 2+ 
ions. 

To be able to account for the antisite disorder in the 
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further analysis of the system, the Cu 2+ /Zn 2+ mixing 
was measured using neutron powder diffraction at the 
Rotax neutron time-of-flight diffractometer at the ISIS 
facility, United Kingdom. We used a 4 g powder sample 
of deuterated x = 1 zinc paratacamite, synthesised using 
the hydrothermal method as described in ll|. The pu- 
rity of the samples and Cu 2+ to Zn 2+ ratio were verified 
with powder x-ray diffraction and Inductively Coupled 
Plasma Auger Electron Spectroscopy (ICP-AES) with an 
accuracy of ±0.03 in x. Neutron-diffraction data were 
collected at 285, 150 and 10 K, and Rietveld analysed 
against the structure as reported in 11]. No structural 
changes were observed with temperature, and the level of 
deuteration was refined to 94.0(6)%. During the refine- 
ment the total Cu 2+ :Zn 2+ ratio was fixed at 3:1, as mea- 
sured directly using ICP-AES. It was further assumed 
there were no vacant Zn or Cu sites. In this way the Cu 2+ 
occupancy on the kagome lattice was refined to 91(2)%, 
corresponding to a Zn 2+ occupancy of the inter-plane 
Zn site of 73(6)% citefootnotel, for the highest-statistics 
data set taken at 10 K. The result is given in figure [TJ 
The refinements of the data taken at 150 and 285 K were 
in overall agreement within the experimental error. On 
relaxation of the constraints the solution was stable, but 
no longer unique. There was a small further reduction of 
the residues for a slight reduction of the Cu and Zn site 
occupancies, which also suggests that the experimental 
error will be larger than stated. 



by Helton et al. [1J| . Figure [2] presents the heat capaci- 
ties of samples with x = 0.8, 0.9 and x — 1 in and 9 T 
and for x = 0.5 in T. For intermediate fields, not shown 
here for clarity, the shoulder gradually moves to higher 
temperatures while the total entropy below ~ 24 K re- 
mains constant. The magnetic susceptibility of samples 
with x — 0.8, x = 0.9 and x — 1 as shown in figure [3] 
was measured with a Quantum Design MPMS system, 
on ~ 50 mg pellets. 
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FIG. 3: The zero- field cooled and field cooled magnetic sus- 
ceptibility for zinc paratacamite with x = 0.8 and x = 0.9, 
compared with the susceptibility for x = 1 which has a negli- 
gible hysteresis. 



□ x = 1,C v (0T)-C v (9T) 
— fit with S=1/2,f = 0.21 
« x = 0.8, C v (0 T) - C v (9 T) 




x=0.5, 0T 
x=0.8, 0T 
o □ x=0.8, 9T 
x=0.9, T 
x=0.9, 9T 
x=1, 0T 
o-ox=1, 9T 



T(K) 



FIG. 2: The heat capacity in T for samples with x = 0.5, 
0.8, 0.9 and 1 as well as in 9 T for x = 0.8, 0.9 and 1. The 
error bars are given for x = 1 only. The inset displays ACv/T 
for x = 0.8 and x = 1 and their respective fits. 



The heat capacity measurements were carried out us- 
ing a Quantum Design PPMS system, on ~ 5 mg 
dye-pressed pellets of Zn 2 ,Cu4_ x (OII) 6 Cl2 with x — 
0.5,0.8,0.9 and 1.0. We could reproduce the heat capac- 
ity for x = 1 in 0, 1, 2, 3, 5, 7 and 9 T fields as reported 



To the eye, the field dependence of the heat capacity 
is similar to the Schottky anomaly arising from defects 
in Zn-doped Y 2 BaNi0 5 and [Ni(C 2 H 8 N2)2(N0 2 )]C104 
(NENP) [2(|. Hence, we have applied a similar analy- 
sis as described in [20(. To study the field-dependent 
part of the heat capacity, for each sample (x) the differ- 
ence was taken between the interpolated heat-capacity 
curves measured in different fields. The inset in fig- 
ure [2] shows the difference between the and 9 T curves 
ACV/T = [C V (H 1 = T) - C V (H 2 = 9 T)]/T for x = 
0.8 (crosses) and for x = 1 (squares). We found that the 
field-dependent part of the heat capacity can be modelled 
by a small number of zero-field split doublets, i.e. inter- 
acting S — 1/2 spins or S = 1/2 excitations. ACV/T was 
fitted with f[Cy =1/2 (AE H1 )-Cy =1/2 (AE H2 )}/T, where 
/ is the fraction of doublets per unit cell (or their spec- 
tral weight). C'y~ 1 ^ 2 (AEh) is the heat capacity from a 
5 = 1/2 spin with a level splitting AEh which for fields 
H > 2 T equals the Zeeman splitting with g w 2.2, as 
shown in the inset of figured] The shoulder in the heat 
capacity in zero-field, which corresponds to a zero-field 
splitting of the doublets of AE ~ 1.7 K (0.15 meV) for 
x = 1, AE ~ 2.1 K for x = 0.9 and AE ~ 2.2 K for 
x = 0.8 indicates that the levels involved are part of an 
interacting system, and cannot be ascribed to a param- 
agnetic impurity phase. The best agreement with exper- 
iment was obtained when a small Gaussian spread a in 
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level splittings AE was taken into account, indicated as 
the error bars in the inset of figure [4] This brings the 
number of fit parameters to 5. The lines through the 
data points in the inset of figure [5] are the fit results for 
x = 0.8 and x = 1. We find that / = 0.21(1), 0.22(1) 
and 0.19(1) for x — 0.8,0.9 and 1.0 respectively. For 
i = l, with three Cu 2+ ions per unit cell, this accounts 
on average for 6.0(6)% of all Cu 2+ . 
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FIG. 4: The field-independent part of the heat capacity as 
obtained from the T and 9 T data of samples with x = 0.8, 
0.9 and 1. The dotted lines (blue online) give fits to the 
corrected heat capacity for x = 1 with jT a with a = 1.3 
and 1.7. The inset shows AE as a function of H between 
and 9 T, compared to the Zeeman splitting with g = 2.2. 
Sample 2 is the deuterated x = 1 sample which was used in 
the neutron diffraction experiment. 



Due to a spin gap the heat capacity of the S = 1/2 
kagome antiferromagnet is expected to show a shoulder 
between J/20 and ~ J/10 corresponding to thepopula- 
tion of the lowest magnetic (S to t = 1) levels (HI. Ift [71 |2~J] . 
In our data a shoulder is evident in zero field. How- 
ever, with the application of a magnetic field this shoul- 
der moves to higher temperatures (energies), as shown in 
the inset of figure 0J This is very different from what can 
be expected for, for example, a singlet-triplet system with 
a non-magnetic ground state. For the latter an applied 
field will lower the energy of the S = 1 level with spin 
aligned along the field, so that for sufficiently strong fields 
a level crossing occurs and this 5=1 level will become 
the ground state. It is clear that such a level crossing 
is not observed here. Most likely, the lowest energy level 
involved in the system giving rise to a field-dependence is 
a magnetic level too. Several models have been tried, of 
which only the doublet (a 5 = 1/2 system) gives an over- 
all consistent fit for all 18 curves from a total of 5 samples, 
for each model using only 5 fitting parameters. A model 
with a triplet of S — 1 levels results in a slightly poorer fit 
to the data as compared with a doublet. Similar models 



with higher- level multiplets (Stot > 1) can not be brought 
into agreement with our data. It should be noted that 
a doublet with a field-dependence as described here has 
also been observed in neutron spectroscopy data 14 , 2^ ■ 
As is shown in figure [31 the system gradually develops a 
magnetic hysteresis as x is lowered (the Cu 2+ concentra- 
tion is increased), while the muon relaxation increase is 
indicative of a slowing-down of the spin dynamics [l3T |. 
The hysteresis is a history dependence which rules out 
a macroscopic quantum state for the system as a whole, 
since such a state would have a unitary time evolution as 
described by the Schrodinger equation. That the latter is 
not the case here is also clear from the energy gap for the 
antisite spins in zero field, which increases as x is low- 
ered. Using our model this increase is quantified as a gap 
of 1.7 K for x = 1 to 2.2 K for x = 0.8. This may be the 
strongest indication that the energy gap corresponds to 
local excitations rather than coherent many-body quan- 
tum states of the total system. In the latter case the 
time dependence should follow the Schrodinger equation 
where a larger gap leads to faster dynamics. 

We suggest, as is also done in 0, E3, EH, that the 
fraction / of zero-field split doublets, which models the 
field-dependence in the heat capacity for 0.8 < x < 1, 
are weakly-coupled 5=1/2 spins from Cu 2+ ions resid- 
ing on inter-plane Zn sites (antisite spins). For x = 1 an 
identical fraction / of Zn 2+ ions must occupy Cu sites on 
the kagome lattice. Once f(x) is known the Cu 2+ cov- 
erage c(x) of the three Cu sites per unit cell is given by 
c = 4 — x — f(x). An important assumption in our argu- 
ment is that the heat capacity of a slightly diamagnet- 
ically doped kagome lattice is field independent, which 
is reasonable as long as g^H <C 9 W 0, [23[. For the 
deuterated x = 1 sample of the neutron-diffraction mea- 
surements it follows that the antisite disorder is 6.3(3)% 
in Cu 2+ or 19.0(9)% in Zn 2+ , in rough agreement with 
the neutron diffraction result. 



TABLE I: The fitted fraction of antisite spins per unit cell 
/, the corresponding Cu 2+ occupancy of the kagome lattice 
c = 4 — x — /, the total entropy from the antisite spins 5/, the 
measured total entropy S(T) up to 24 K and the percentage 
of the entropy recovered per Cu 2+ spin. 
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±0.03 


/ 

±0.018 


c 

±0.03 


S f /R 
/ln(2) 


S{T) /R 
<ST = 24 K 6 


S(T) 
(4,-x) ln(2) 
/% 


0.50 


0.50 a 


3.00 


0.346(8) 


1.061(12) 


43.9(2) 


0.80 


0.210 


2.97 


0.15 


0.993(11) 


44.7(2) 


0.90 


0.220 


2.88 


0.15 


0.959(9) 


44.8(2) 


1.00 


0.190 


2.81 


0.13 


0.933(9) 


44.8(2) 



a Here / was not obtained from the heat capacity which at this 
level of doping is altered due to a cooperative transition. 
However, for x = 0.5 it is safe to assume that c = 3.0 (full 
occupancy) and hence f = x. b For T > 24 K no relative changes 
in <S(T) occur between samples with different x. 

Comparing the heat-capacity data of several x = 1 
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samples, all synthesised at a temperature of 484 K, an 
average antisite disorder of ~ 6.0(6)% in Cu 2+ is derived, 
as listed in table U along with the results for x = 0.8 and 
x = 0.9. The chemical potential behind the Cu 2+ /Zn 2+ 
partitioning can now be estimated to ~ 1400 K, a plau- 
sible value given that most likely the Zn site becomes 
locally slightly angle-distorted, if occupied by an other- 
wise orbitally degenerate Cu 2+ ion. The Cu sites on the 
kagome lattice are energetically favoured by the Cu 2+ 
ions, and there is only a slow increase of the Cu 2+ oc- 
cupancy on the Zn sites until the Cu 2+ occupancy of 
the kagome lattice (c in table [TJ) is almost complete. For 
x = 0.8 the magnetic hysteresis is too large to be as- 
cribed to impurities or local variations in Zn stoichiom- 
etry (Fig. [3]). Since even at x = 0.8 only ~ 20% of the 
Zn sites are occupied by Cu 2+ , this hysteresis must be 
due to the higher connectivity of the 3D lattice, which is 
mainly due to the higher Cu 2+ occupancy of the kagome 
planes (see table HJ. Hence, for the phases x < 1 which 
have a magnetic hysteresis the kagome layers must be in 
a magnetic state, i.e. both singlet and triplet states mix 
into the ground state. Since no quantum phase transi- 
tion occurs between x = 0.8 and x — 1 as is clear from 
our heat capacity data, the ground state of the kagome 
layers in the x = 1 phase must be mag netic too. This is 
in support of NMR measurements 2J, [25| in that there 
is no spin gap. What is remarkable in the present case 
is that the appearance of unpaired spins precedes the 
breaking of spin-rotational symmetry to a long-range or- 
dered state. This results in a macroscopically degenerate 
ground state with increasingly glassy dynamics as x is 
lowered. 

The heat capacity of the kagome lattice can be esti- 
mated by subtracting the heat capacity from the Cu 2+ 
spins on the Zn sites. The result for the data with 
0.8 < x < 1 is shown in figure [4] For all x the curves 
obtained from the and 9 T data are identical within 
the experimental error, which follows from the quality 
of the fit as described in the previous paragraphs. This 
part of the heat capacity most likely corresponds to the 
kagome layers. In this field-independent part of the heat 
capacity a weak shoulder is visible at a slightly higher 
temperature than the shoulder due to the antisite spins. 
As x — > 1 the shoulder becomes less pronounced, and 
hence, it may be interpreted as due to the entropy re- 
lease when the fluctuations in neighbouring kagome lay- 
ers, which are connected via the antisite spins, decouple. 
If the heat capacity from perfectly 2D kagome layers fol- 
lows a power-law for T — > 0, then taking into account 
the distortive effect of the shoulder due to couplings be- 
tween the kagome layers, O.IT" mol -1 formula units with 
a = 1.3(1) is our best estimate. An exponent a = 2 can- 
not be brought into agreement with our data. 

In summary, based on the field-dependence of the 
shoulder in the low-temperature heat capacity, which cor- 
responds to the weakly dispersive feature observed at the 



Zeeman energy in neutron spectroscopy data, we rule 
out interpretations based on a singlet-triplet splitting. 
Rather, the feature remains a doublet over the entire 
range of applied fields. We suggest these doublets are the 
magnetic states of the Cu 2+ antisite spins, which raises 
an important question as to the origin of the observed 
zero-field splitting of the antisite spins. From analysis 
of the heat capacity and the magnetic susceptibility as 
a function of x we further conclude that even for x = 1 
the ground state of the kagome system is a gapless spin 
liquid. 
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